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ABSTRACT 

BL Lac objects of the intermediate subclass (IBLs) are known to emit a 
substantial fraction of their power in the energy range 0.1-10 GeV. Detecting 7- 
ray emission from such sources provides therefore a direct probe of the emission 
mechanisms and of the underlying powerhouse. 

The AGILE 7-ray satellite detected the remarkable IBL S5 0716+714 (z ~ 
0.3) during a high state in the period from 2007 September - October, marked by 
two very intense flares reaching peak fluxes of 200 x 10~ 8 photons cm _2 s _1 above 
100 MeV, with simultaneous optical and X-ray observations. We present here a 
theoretical model for the two major flares and discuss the overall energetics of 
the source. 

We conclude that 0716+714 is among the brightest BL Lac's ever detected at 
7-ray energies. Because of its high power and lack of signs for ongoing accretion 
or surrounding gas, the source is an ideal candidate to test the maximal power ex- 
tractable from a rotating supermassive black hole via the pure Blandford-Znajek 
(BZ) mechanism. We find that during the 2007 7-ray flares 0716+714 approached 
or just exceeded the upper limit set by BZ for a black hole of mass 10 9 M . 

Subject headings: BL Lacertae objects: individual (S5 0716+714) - gamma rays: 
observations 
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1. Introduction 

Blazars constitute a class of active galactic nuclei (AGN) that often show very strong 
and rapid flux variability over the electromagnetic spectrum. They are widely held to 
contain a black hole (BH) with mass in the range 10 7 — 1O 9 M , that launches relativistic 
jets emitting highly non-thermal radiation. 

The jet transports energy in electromagnetic form, and bulk plus random kinetic energy 
of charged particles. The source radiation may also show a contribution by the accretion 
disk (including the big blue bump, BBB), the broad line region (BLR), and a dusty torus 
(see Urry & Padovani 1995); lack or weakness of such contributions mark out the class of 
BL Lac objects. The primary energy of the jet may be supplied by power extracted from 
the central rotating BH via interaction with its accretion disk (Blandford & Znajek 1977; 
Blandford & Payne 1982; see also discussions in Cavaliere & D'Elia 2002 and McKinney 
2005). 

The blazar 0716+714 is a distant BL Lac at z = 0.31 ± 0.08 (Nilsson et al. 2008); 
its optical-UV continuum is so featureless (Biermann et al. 1981; Stickel, Fried & Kiihr 
1993) that a redshift estimate has been possible only resolving and using the host galaxy 
as a standard candle (Nilsson et al. 2008). This BL Lac is of the IBL type, displaying a 
first broad peak in the optical-UV bands and showing another broad peak near 1 GeV; the 
crossover between the two components falls in X-ray range (Massaro et al. 2008; Ferrero et 
al. 2006). Recently, 0716+714 has been detected by AGILE above 100 MeV several times 
during the period September-October 2007 when the source was quite active in optical 
band with variations on 1-day timescale (Villata et al. 2008). Two bright 7-ray flares 
were detected. The first one reached a flux of ~ (200 ±40) x 10 _8 photonscm _2 s _1 on 
2007 September 11 with a photon index of ~ 1.6 and duration < 1 day. In 2007 October, 
following another prominent optical activity detected by the WEBT consortium, AGILE 
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and Swift satellites pointed again at the source. On 2007 October 23, another 1-day bright 
flare was observed with a 7-ray flux comparable to that detected in September (Chen et al. 
2008). Around this date a bright one-day flare and strong day- variability were observed 
in the optical band (Villata et al. 2008); UV and soft X-rays also showed strong and fast 
variability, whereas modest or none variability appeared in the band 4-10 keV (Giommi 
et al. 2008). Meanwhile, the radio flux showed a slow coherent increase that remarkably 
begins around the day of the first 7-ray flare and culminates around the date of the second 
7-ray flare (Villata et al. 2008). Thereafter 0716+714 was detected going back to its ground 
state in all bands, with the 7-ray photon index that softened toward 1.9 (Chen et al. 2008) 
as previously observed by EGRET (Lin et al. 1995). 

In this paper, we present a physical model for the 0716+714 flaring states, and study 
their extreme energetics with implications on energy extraction from a rotating BH. 



2. Spectral modeling 

2.1. One Simple Model, Synchrotron Self-Compton 

The simplest homogeneous synchrotron self-Compton (SSC) model assumes the blazar 
emissions to be produced in a "blob" of radius R, containing relativistic electrons in a 
combination of tangled and uniform magnetic field. The emitters move toward the observer 
with bulk Lorentz factor T (see, e.g., Tavecchio, Maraschi & Ghisellini 1998). We assume 
the emitters to emerge from the injection/acceleration phase with a jet-frame distribution 
of the random energies 7772c 2 in the form of a standard broken power-law 

Ue{l) = (7/7^ + (7/7^ 2 ' (1) 
where Ci and ( 2 are the spectral indices for 7 < % and 7 > %, respectively, % is the 
Lorentz factor at the break. These electrons emit a primary synchrotron spectrum; a second 
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contribution is then produced by inverse Compton (IC) as the primary synchrotron photons 
scatter off the same electron population. The spectral energy distribution (SED) behaves 
as eF(e) oc e 1_a , where e is the energy of the received photons, and a = (( — l)/2. 

For electrons in a magnetic field B, the synchrotron SED peaks around 

where h is Planck's constant, z is the redshift of the source, and 5 = [T(l — (3 cosO)]' 1 is the 
bulk Doppler factor due to the flow of emitters toward the observer at an angle 9 relative 
to the line of sight; the SED at the synchrotron peak is 

e s F(e s )cxS 4 R 3 B 2 K^ 2 . (3) 

As to the IC component, its SED contribution peaks at 

<« = ^ (4) 

with a peak value of 

e c F(e c )oc5 4 J R 4 J B 2 K 2 7 fe 4 (5) 

if the scattering takes place in the Thomson regime with the density of target photons 
scaling as n p h oc F s R/c. The relativistic motion toward the observer amplifies the emitted 
power by the factor <5 4 , and allows it to vary on a timescale 

Ml + z) 

Uar > ^ W 

close to or shorter than the crossing-time t cr = R/c. 

Due to the synchrotron and IC losses the electrons cool with timescale 

3mc 



Tcoold) 



4P 2 <TTl(U B + U r ) 

where ot is the Thomson cross section, U B — B 2 /8n and U r is the energy density of radiation 
before scattering. This sets a typical cooling break at 7 C00 2 = 3mc 2 /4<7t R(3 2 (Ub + U r ) 
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beyond which the electrons cool rapidly. In the following, we take into account this 
constraint on the particle distributions. 

Simultaneous multi-frequency observations can provide the five quantities R, 5, B, K e 
and 7b with the five Equations (2)- (6). 

For BL Lacs with high frequency peaks (HBLs) requiring electrons of higher 
energies (75 > 10 4 ) the scattering approaches the Klein-Nishina (KN) regime with a 
blob-frame photon energy > m e c 2 /7&. In the extreme KN regime the IC SED peaks at 
e c ~ 7i>m e c 2 <5/(l + z); the dependence on B and 75 progressively weakens as the two latter 
parameters grow. 



2.2. An Addition to the Model, External Seed Photons 

Additional target photon can be provided by a source external to the jet (see Dermer 
et al. 2009). In this case, the high-energy component of the spectra is due to the electrons 
that Compton-scatter the external photons (EC); the SED now peaks at energies 

and the corresponding SED value is 

e c F(e c )^5'R 3 K 1 2 b N' ext e' ext . (8) 

In this EC process two new ingredients enter: e' ext and N' ext , respectively, the energy 
and the density at peak of the external photons as seen by the moving blob. This has two 
main consequences: 

i) The model contains two further degrees of freedom and the parameter evaluation may be 
degenerate; 

ii) These news quantities are related to N ext and e ext in the observer frame by means of the 
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bulk Lorentz factor T in a manner that depends on the geometry of the system (Dermer & 
Schlickeiser 2002), causing an additional dependence on T in the EC spectra. Dermer & 
Schlickeiser (1993) discuss SED dependences on T varying from oc T 3 to oc T 6 , for photons 
entering into the blob from behind or head-on, respectively. 

2.3. Flux Variation Patterns 

Equations. (2) - (8) show that, in the synchrotron-IC framework, variabilities of the 
first and second peaks are correlated, possibly with a lag t de i ~ t cr (l + z)/S. 

When the energy fluxes around these peaks oc eF(e) are simultaneously monitored, 
we can compare the corresponding light curves <p{t) in the respective band energy 
e s and e c , possibly with a time-lag tdei- Then given two times t\ and £2, the ratio 
r c = 4> c {t2 + tdei) / <p c (ti + tdei) shows a dependence on r s = (frsfo) / <f>s(ti) that is related to 
the emission mechanism. Here we consider some relevant cases. 

If the variability is mainly due to electron injection/acceleration, K and/or 7^ changes; 
then, r c = r 2 s results if SSC dominates the IC emission. If instead EC dominates then r c = r s 
applies (compare Equation (3) with Equations (5) and (8)). If variations are mainly due to 
changes in B, then r c = r s holds for SSC, and r c — 1 applies for EC. If instead T varies, we 
have r c = r s in the SSC; in the EC framework we have different behaviors depending on the 
geometric of the jet relative to the external radiation: r c = r x s with x ranging from 3/4 for 
seed photons entering from behind, to 3/2 for photons entering head-on the flow. 

In Table 1 we report the relations for the Thomson and the extreme KN regimes (see 
also Paggi et al. 2009). The cooling of electrons acts as a variation of K and 7^ Hence, in 
a flare due to electron injection/acceleration, the trajectories r c versus r s remain unchanged 
by pure radiative cooling. 
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2.4. Modeling the 7-ray flaring states 

A 1 day time lag between the emission in optical and 7-ray bands during the two flares 
is apparent from Figures 1 and 3 of Chen et al. (2008). This constrains the emitting region 
radius and the Doppler factor to R < 5 x 10 16 (5/20) cm; the duration of both flares is 1 
day or less, which argues for cooling times r coo ;(7t,) ~ R/c. 

Moreover, optical and gamma light curves around the two flare dates show evidence 
that r c = r 2 s applies (see Figure 2 and its caption): this argues for a SSC process in the 
Thomson regime, and concurs with the lack of sign of external gas to rule out EC process 
(see also Table 1). 

Previous radio monitoring by the Very Long Baseline Array (VLBA) telescope of 
0716+714 showed the presence of more super-luminal components relative to the active 
state of 2003-2004 (Rastorgueva et al. 2009). Moreover, the absence of the signatures of 
IC catastrophe provided a lower limit 5 > 14 for the Doppler factor (Ostorero et al. 2006; 
Fuhrmann et al. 2008; see also Wagner et al. 1996); Bach et al. (2005) also argue for high 
Doppler factors in these fast variable components and for a viewing angle 9 < 5°. 

We note that modeling the SED of 0716+714 with a standard one-zone SSC model 
would fail to reproduce the simultaneous radio, optical, X-ray and 7-ray data for the 
October flare. In particular, the crossover in the X-ray band would not be well reproduced, 



Table 1: Variability Patterns. 



Model 


r changes 


B changes 


K changes 


7b changes 


SSC Th. 
EC Th. 


r c = r s 
r c = r s x 


r c = r s 
r c =l 


r c = r s 2 
r c = r s 


r c = r 2 
r c = r s 


SSC KN 
EC KN 


r c = r s 
r c = r s x 


r c -» r s °- 5 
r c =l 


r c = r s 2 
r c = r s 


r c -> 1 
r c -> 1 
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and the hard X-ray flux would be overestimated by a factor ~3; furthermore it would be 
difficult to model the hardness of the 7-ray spectrum during the September flare (Figure 1 
red dashed lines). A one-component model also hardly explain together the slow trends of 
the radio, optical and hard X-ray bands and the faster variability observed in the optical, 
soft X-ray and 7-ray bands (see Villata et al. 2008, Giommi et al. 2008, Chen et al. 2008). 
Hence, we adopt a two-component model: the first produces the slowly variable radio and 
hard X-ray bands, whereas the second is responsible for the faster variability in optical-UV, 
soft X-, and 7-ray bands. Despite these problems, the possibility to fully constrain by 
observations the parameters stimulate us to also show for comparison a one-component 



The SED of all these models are shown in Figure 1 and the parameters are listed in 
Table 2; a viewing angle 6 2° is adopted according to Bach et al. (2005). 



Under the assumption of isotropic emission, the observed power radiated from a source 
with luminosity distance Dl(z) is 



The jet transports a total power Ptot, flare = L r + L kin + L B contributed by intrinsic 
radiated power, kinetic energy flow of the electrons and of the cold protons (with one proton 
per emitting electron), and Poynting flux, respectively: 



model. 



3. The extreme energetics of 0716+714 




L obs T 2 /5 A = 2 x 10' 3 (5/15)- A T 1 2 L obs , 



(9) 



L e = 



0.8 x 10 43 Rl 



i 6 Ti 2 (n e7 ) 4 ergs 



(10) 
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L p = 1.4 x lO^i&Ti 2 {n p ) 1 ergs" 1 , (11) 
L B = 3.8 x lO 43 ^ Ti 2 B 2 ergs" 1 , (12) 

see also Celotti & Ghisellini (2008). The latter authors show that in BL Lacs L r tends to 
match the sum of the other contributions. In fact, for the two flares of 0716+714 at redshift 
z = 0.31 we find L r ~ 2 x 10 45 erg s -1 , and from our two-component SSC model (with the 
parameters listed in Table 2) we obtain a total jet power 

Pto^are = (3.5 ±l)xl0 45 ergs" 1 

with L r > (L e + L p + L B ). Under this condition the total jet power is minimized and 
the details of cooling does not affect materially the global energetics, being the radiated 
luminosity L r mainly contributed by peaks emission. Moreover, the uncertainty in Ptot, flare 
is mainly due to the observed 7-ray flux error. 

For the one-component model, we obtain Ptot, flare — (1 ± 0.5) x 10 46 ergs _1 and 
L r < (L e + L p + Lb) holds. In this case, the total jet power is not dominated by the radiated 
one, but the parameters now are well constrained by the observation and the uncertainty is 
still due to the flux error. 



3.1. Testing the Blandford-Znajek Mechanism 

Here, we compare the jet powers provided by our models with the BZ mechanism; this 
set a limit on the power extractable from a rotating BH 

P BZ ~ 2 x lO^Mgergs" 1 (13) 

under conservative values of B, as discussed in Cavaliere & D'Elia (2002, and references 
therein; see also our discussion). 
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Estimating the redshift and the BH mass of 0716+714 is not trivial because of the lack 
of emission lines. Recently, however, Nilsson et al. (2008) pinpointed the host galaxy of 
this source and reported a determination of its redshift at z = 0.31 ± 0.08. Then, for this 
host, a BH of mass M ~ 5 x 1O 8 M should accord with the fundamental plane of BL Lacs 
(Falomo et al. 2003). 

On the other hand, some authors using micro-variability of the optical flux estimate 
the mass of the central BH for 0716+714 by 

M < °] Tb r « 3 x 1O 8 M [— — 1. (14) 
6G{l + z) OL 20 600 J V ; 

Considering r ~ 450 s as the shortest variability time found by Sasada et al. (2008) one 

obtains M < 2 x 1O 8 (5/2O)M . Gupta, Srivastava & Wiita (2009) obtain the more stringent 

constraint M < 5 x 1O 7 M for the case of a Kerr BH. However, these authors discuss that 

the micro- variability in 0716+714 may be due to a small region in the jet or due to internal 

disk modes, so that the BH mass would be higher. 

We note that a BH of mass M < 1O 8 M in 0716+714 would imply in Equation (13) a 
power limit Pbz « L r inconsistent with the power L r emitted during the flares. 



4. Discussion and Conclusions 

We modeled the 0716+714 flares of 2007 September 11 and October 23 with a 
two-component SSC model, similarly to Tavecchio & Ghisellini (2009). Despite the larger 
number of model parameters, we believe that our two-component modeling reproduces the 
complex variability and the hard 7-ray spectra of 0716+714 better than a one-component 
model. 

Our Figure 2 indicates a quadratic dependence between the synchrotron and IC fluxes. 
This concurs with the lack of emission lines and BBB to rule out models involving external 



_o -in 1 1 1 1 1 1 1 1 1 1 1 1 r 

10 12 14 16 18 20 22 24 26 




log [V / Hz] 



Fig. 1. — Bottom panel: SED relative to the 2007 September 11 data. Curves labeled "c 
I" and "c II" (gray and blue solid lines) show the two separate components. Data in light 
gray represent a previous low state plus the EGRET observation in 7-ray (see Lin et al. 
1995). Top panel: two- components SED relative to the 2007 October flare (solid line), and 
after 1 day by radiative cooling (dotted line). In both panels black stars are simultaneous 
data (see Villata et al. 2008; Giommi et al. 2008; Chen et al. 2008), the black lines are the 
two- components models and the dashed red lines represent the best one-component models. 



Fig. 2. — Variability plane with the trajectories of 0716+714 based on simultaneous optical 
and 7-ray observations. Open circles are for 2007 September, and filled triangles for 2007 
October data. The thick solid line represents the quadratic trend, expected from electrons 
injected or accelerating in the SSC framework; for comparison, the linear behavior expected 
from EC is represented by the thin dotted line. 
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sources of seed photons that produce a linear dependence, and to strongly support electrons 
radiating in the Thomson regime within the SSC framework as the most likely radiation 
process. We refer to our Figures 2 and 1 in Chen et al. (2008) to show that the 2007 
September data lie around the maximum; in time, the trajectory starts from the lower r c 
value, attains a maximum and then falls down. In 2007 October, the trajectory describes 
the decline of a flare starting from the higher r c value. Rises and falls both occur on 1-day 
timescale, and the trajectories are quadratic. This suggests the electron acceleration and 
cooling to occur on similar timescales. 

The crossover between the synchrotron and IC branches of the SED is located in X-rays 
(see also Foschini et al. 2006; Ferrero et al. 2006), and is contributed by both the slowly 
variable component (I) and the faster component (II). 

As to the faster components (II) adopted for the two flares, they are marked by high 
electron energies % ~ 7 x 10 3 with a sharp low energy cutoff 7 mi „ ~ 2 10 3 (see also Tsang 
& Kirk 2007). Moreover, high bulk Lorentz factor T = 15 (that is 5 ~ 23) is used in 
accord with Wagner et al. (1996). In Table 2, it is shown that our choice of parameters 
implies r coo i{j h ) ~ t cr for the fast components II: those quench very rapidly causing strong 
variability in the optical-UV, soft X- and 7-ray bands. Little or no variability results in 
radio and hard-X ray bands produced by the rising part of the slow components I. This 
behavior is in agreement with the complex multi-band variability reported by Chen et al. 
(2008), Giommi et al. (2008), and Villata et al. (2008). 

Considering the redshift z = 0.31 of 0716+714 , we find that its intrinsic radiative 
luminosity is of order L r w 2 x 10 45 erg s _1 . On adding the other jet components in 
Equations (10) - (12), the total power becomes Ptotjiare ~ 4 x 10 45 erg s" 1 for the 
two-component model, and P to t, flare ~ 2 x 10 46 erg s" 1 for the one-component model: the 
source exceeds the BZ limit Pbz ~ 2 x 10 45 M B h/M 9 erg s" 1 (see Figure 3). This obtains 
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Fig. 3. — Observed high-energy peak fluxes (top panel) and the corresponding intrinsic 
peak luminosities (bottom panel) for a number of BL Lac objects during their flaring states, 
ordered with their z: Mrk 421 (open symbols, here shown for reference), BL Lac, W Comae, 
and 0716+714 . For the first three objects, data are respectively from Donnarumma et al. 
(2009), Ravasio et al. (2002), and Bottcher et al. (2002). The shaded area represents the 
BZ limiting luminosity range for a BH mass in the range (3 x 10 8 — 10 9 ) M & . The arrow 
points at the total jet power by also adding the components of Equations (10) - (12). 
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for a maximally rotating BH of lO 9 M (spun by past accretion) via interaction with a 
disk magnetic field B ~ 10 4 G sustained by gas or radiation pressure, in the absence of 
ongoing accretion (see discussion by Cavaliere & D'Elia 2002). The simultaneous nature of 
our multifrequency data during the 2007 7-ray flares of 0716+714 is a crucial ingredient 
of our result. Other BL Lacs with weak or negligible accretion disk or BLR contributions 
have been reported to attain (model dependent) high total luminosities, e.g., 2032+107 
with an inferred P ~ lO^ergs" 1 (see Celotti & Ghisellini 2008). Corbel & Reyes (2008, 
ATel 1744) report the high 7-ray peak flux recently attained by 0235+164 at z ~ 0.94 for 
which, however, EC contributions cannot be ruled out. However, to our knowledge, only 
the simultaneous data obtained by our group for 0716+714 provide a model independent 
evaluation of the total jet power from a BL Lac being close or just above the BZ limit. 

We also show in Figure 3 the intrinsic radiated power for other BL Lacs with intense 
7-ray and TeV emissions as for the 1997 flare of BL Lacertae with peak flux around 
500 x 10~ 8 photons cm~ 2 s~ 1 but with lower redshift z = 0.069 (Bloom et al. 1997; see also 
Ravasio et al. 2002) and intermittent evidence of lines and thermal emissions. We also 
report W Comae with redshift z = 0.102 but 7-ray flux at levels of 50 x 10~ 8 photons cm _2 s _1 
(Bottcher, Mukherjee & Reimer 2002). 

The increasing trend of L r with z is likely to arise from the Malmquist bias, while the 
decrease at low z may result from sampling a limited cosmological volume. Nevertheless, 
we stress that up to now no BL Lac source has sharply exceeded the BZ limiting power. It 
will be worthwhile to keep under close watch the most distant BL Lacs (despite the obvious 
monitoring difficulties) to catch powerful emissions. If violations will be found, these may 
be possibly discussed in terms of the Blandford-Payne mechanism (Blandford & Payne 
1982) that, however, requires ongoing accretion not supported in the case of 0716+714 . 
Alternatively, higher powers may be attained with higher magnetic fields up to B 2 /An < pc 2 
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related to the plunging orbits (see Meier 1999), which however imply short source lifetimes 
in the absence of accretion. 

In conclusion we find the total power transported in the jet of 0716+714 to be 
Ptot, flare ~ 4 x 10 45 ergs _1 , and so to approach the limit of the BZ mechanism for a BH 
up to 10 9 M Q with conservative values of B. Such high powers in 0716+714 constitute an 
unescapable consequence of two observed facts: 

• The 7-ray flux attains high levels in the SED during the two flares, with an emitted 
power comparable with the optical one as shown by simultaneous observations (see 
Figure 1); 

• The lack of external sources of seed photons related to emission lines or BBB concurs 
with the observed quadratic dependence r c = r\ to rule out EC contributions to the 
high energy hump (see Figure 2) and considerable ongoing accretion. 

Acknowledgments: This investigation was carried out with partial support under ASI 
contract no. 1/089/06/2. 
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Table 2: Model Parameters for the 2007 7-ray Flares of 0716+714 . 



Date 


Comp. 


r 


B (Gauss) 


R(cm) 


K(cm 3 ) 


7b 


7min 


Ci 


C2 


T CO ol(7b)/t C r 


Sept 11 


c I 


10 


0.4 


3 x 10 16 


3.5 


3800 


100 


2 


4.5 


1.3 




c II 


15 


0.3 


3.5 x 10 16 


1.4 


7000 


3500 


2 


5 


1.0 




single 


15 


0.3 


3.5 x 10 16 


1.8 


7000 


60 


1.8 


5 


1.0 


Oct 23 


c I 


10 


0.4 


3 x 10 16 


2.5 


4000 


40 


2 


4.5 


1.2 




c II 


15 


0.3 


3.5 x 10 16 


1.5 


6500 


1300 


2 


5 


1.1 




single 


15 


0.3 


3.5 x 10 16 


1.8 


6500 


30 


1.8 


5 


1.1 
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